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The orientation of galaxy dark matter haloes around cosmic voids 
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ABSTRACT 

Using the Millennium N-body simulation we explore how the shape and angular momen- 
tum of galaxy dark matter haloes surrounding the largest cosmological voids are oriented. We 
find that the major and intermediate axes of the haloes tend to lie parallel to the surface of the 
voids, whereas the minor axis points preferentially in the radial direction. We have quantified 
the strength of these alignments at different radial distances from the void centres. The effect 
of these orientations is still detected at distances as large as 2.2 R VO id from the void centre. 
Taking a subsample of haloes expected to contain disc-dominated galaxies at their centres 
we detect, at the 99.9% confidence level, a signal that the angular momentum of those haloes 
tends to lie parallel to the surface of the voids. Contrary to the alignments of the inertia axes, 
this signal is only detected in shells at the void surface (1<R<1.07 R VO id) and disappears 
at larger distances. This signal, together with the similar alignment observed using real spi- 
ral galaxies (Trujillo, Carretero & Patiri 2006), strongly supports the prediction of the Tidal 
Torque theory that both dark matter haloes and baryonic matter have acquired, conjointly, 
their angular momentum before the moment of turnaround. 
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1 INTRODUCTION 

The angular momentum of a galaxy plays a central role in deter- 
mining its evolution and final type. However, understanding the 
origin and properties of the galactic angular momentum has been 
one of the key problems in astrophysics in the last five decades. 
The current 'standard' theory for the origin of the angular momen- 
tum, within the cosmological framework of hierarchical structure 
formation, is the tidal-torque theory (hereafter TTT). This theory 
is based on early ideas from Hoyle (1951) that suggested that the 
angular momentum of a galaxy arises from the tidal field of neigh- 
bouring galaxies. This idea was further developed and quantified 
by Peebles (1969), Doroshkevich (1970) and White (1984). 

The TTT suggests that most of the angular momentum is 
gained gradually by the protohaloes during the linear regime of 
the growth of density fluctuations, due to tidal torques from neigh- 
bouring fluctuations. Angular momentum grows linearly with time 
at this early epoch and saturates when the halo decouples from 
the expanding background at the moment of turnaround. It is as- 
sumed that during this phase the baryonic component follows the 
dark matter distribution and consequently gains a similar specific 
angular momentum to that of the halo. A subsequent large collapse 
factor of the baryonic matter to the centre of the halo will explain 
the centrifugally supported nature of the galactic discs (Fall & Efs- 
tathiou 1980). 

To first order, the angular momentum of the haloes results 



from the misalignment between the principal axes of the inertia mo- 
mentum tensor (Iy) of the matter being torqued and the principal 
axes of the shear or tidal tensor {Tij=-d 2 qb/dxidxj) generated by 
neighbouring density fluctuations. The leading term of the torque 
is given by LiOcTjfc^jj-Ifcfc), where i, j, and k are cyclic permuta- 
tions of 1, 2 and 3. Several aspects of this picture have been con- 
firmed by a number of studies of the angular momentum properties 
of dark matter haloes both analytically (Heavens & Peacock 1988; 
Catelan & Theuns 1996) and in N-body simulations (Barnes & Ef- 
stathiou 1987; Sugerman, Summers & Kamionkowski 2000; Lee & 
Pen 2000; Porciani, Dekel & Hoffman 2002a,b). 

If I and T were uncorrelated (as frequently has been assumed 
under the argument that the former is local and the latter is domi- 
nated by the distribution of the large-scale structure) the direction 
of the angular momentum in the linear regime should be aligned 
with the intermediate axis of I (the direction that maximizes the 
difference between Ijj and Ikk)- However, Porciani et al. (2002b) 
have found in their simulations that there is a strong correlation 
between the I and T tensors, in the sense that their minor, major 
and intermediate principal axes tend to be aligned at the protohalo 
stage. The strong correlation between I and T should produce an 
angular momentum vector of the haloes that is perpendicular to the 
minor axis of the sheet they are embedded in (i.e. perpendicular to 
the direction of the maximum compression of the large-scale struc- 
ture at that point). However, this last correlation, at least at redshift 
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zero, is expected to be very weak (or totally erased) by non-linear 
effects (such as exchange of angular momentum between haloes) 
at late times. 

From the observational point of view, Trujillo, Carretero & 
Patiri (2006) have found that the rotation axes of spiral galaxies lo- 
cated on the shells of the largest cosmic voids lie preferentially on 
the void surface (in agreement with the expectation for the angu- 
lar momentum of haloes given in the previous picture). The obser- 
vational relation could be explained then as a consequence of the 
spin of the baryonic matter still retaining memory of the angular 
momentum properties of the haloes at the moment of turnaround 
(Navarro, Abadi & Steinmetz 2004). It is key, consequently, to ex- 
plore whether the signal is also found in the haloes of cosmological 
N-body simulations when we mimic the observational technique. 
If so, this will strongly support our current understanding of how 
haloes and baryonic matter have acquired, conjointly, their angular 
momentum. 

Together with the orientation of the angular momentum, 
the alignment of the shape of the galaxy dark matter haloes 
(M<10 13 /i _1 Mq) with their surrounding large-scale structure can 
have important observational consequences. In fact, the shapes of 
dark matter haloes can affect the coherence of tidal streams (Sack- 
ett 1999), can be related to galactic warps (Ostriker & Binney 1989; 
Debattista & Sellwood 1999; Lopez-Corredoira et al. 2002) or can 
affect the distribution of the orbits of satellite galaxies (Holmberg 
1969; Zaritsky et al. 1997; Sales & Lambas 2004; Agustsson & 
Brainerd 2006; Yang et al. 2006). In addition, infall of material 
into the haloes is not isotropic but is expected to be through the 
filaments where the haloes are embedded. Consequently, the orien- 
tation of the dark matter haloes within these structures can affect 
the characteristics of the galaxy properties previously mentioned. 

The alignment of massive (group and cluster) haloes 
(M> 10 13 /i~ 1 Mq) with their surrounding large-scale structure has 
been explored in detail (Splinter et al. 1997; Onuora & Thomas 
2000; Faltenbacher et al. 2002; Kasun & Evrard 2005; Hopkins et 
al. 2005; Basilakos et al. 2006). These works indicate that the ma- 
jor axes of neighbouring galaxy clusters are aligned. A result that is 
in agreement with the "Binggeli (1982) effect". The origin of these 
alignments is still under debate and could be associated to infall of 
material (van Haarlem & van de Weygaert 1993) and/or tidal fields 
(Bond et al. 1996). Due to the lack of resolution in previous sim- 
ulations, it is only now that an exploration of the alignment of the 
galaxy dark matter haloes has become possible (Bailin & Stein- 
metz 2005, hereafter BS05; Altay, Colberg & Croft 2006, here- 
after ACC06; Patiri et al. 2006b, hereafter PA06; Aragon-Calvo 
et al. 2006). These works suggest that the major axis of the haloes 
lies along the filaments. The quantification of the strength of these 
alignments is key to studies of strong and weak lensing where the 
intrinsic distribution and alignment of galaxy shapes plays an im- 
portant role in interpreting the signal (see e.g. Heavens, Refregier 
& Heymans 2000; Croft & Metzler 2000; Heymans et al. 2006). 

The aim of this paper is to characterise the alignment of both 
the shape and angular momentum of galaxy dark matter haloes with 
their surrounding large-scale structure to an unprecedent statistical 
level using the Millennium simulation (Springel et al. 2005). In par- 
ticular, we will focus our attention on haloes surrounding the largest 
cosmological voids. Contrary to filaments (which are strongly af- 
fected by redshift-space distortion), large cosmological voids are 
a feature easy to characterise from the observational point of view 
(Trujillo et al. 2006). In addition, another important advantage of 
the void scheme is that (because of the radial growth of the voids) 
the vector joining the centre of the void with the galaxy (halo) is a 



good approximation to the direction of the maximum compression 
of the large-scale structure at that point. 

This naturally generates a framework for exploring the align- 
ments of the shape and angular momentum of the haloes with their 
surrounding matter distribution. Consequently, our work mimics 
the observational framework to provide an easy interpretation of 
current and future observations. The large volume sampled by the 
Millennium simulation combined with the excellent spatial resolu- 
tion allows us to conduct this project. This unprecendent statistical 
power is absolutely crucial if we want to explore signals expected 
to be very weak like the alignment of halo angular momentum and 
the large-scale structure. 

This paper is structured as follows. Section 2 provides a de- 
scription of the Millennium simulation itself, the void and halo 
samples used, and details of how the shapes and spins were mea- 
sured. Section 3 describes our results and we discuss our findings 
in Section 4. 



2 NUMERICAL SIMULATION 
2.1 N-body simulation 

The main simulation we have used for this study is the Millen- 
nium Simulation of Springel et al (2005). This employs 2160 a 
dark matter particles each of mass 8.6 x 10 8 /i _1 Mq within a 
comoving box of side 500/i _1 Mpc. This simulation was per- 
formed in a ACDM universe with cosmological parameters: Qa = 
0.75, Q, M = 0.25, D. b = 0.045, h = 0.73, n = 1, and a 8 = 0.9, 
where the Hubble constant is characterised as 100/ikms _1 Mpc _1 . 
These cosmological parameters are consistent with recent com- 
bined analyses from WMAP data (Spergel et al 2003) and the 2dF 
galaxy redshift survey (Colless et al 2001), although the value for 
as is a little higher than would perhaps have been desirable in ret- 
rospect. The spatial resolution is 5h _1 kpc everywhere inside the 
simulation volume. 



2.2 Void and halo samples 

To explore the effect of the large-scale structure on the orienta- 
tion of the dark matter haloes, we use a similar technique to the 
one adopted in analysing the orientation of disc galaxies in real ob- 
servations (Trujillo et al. 2006). To this aim we need, first, to find 
and characterise the radius of the large voids in the simulation. To 
speed up the process of finding the voids, we used as an initial guess 
for the position of the void centres the position of the most under- 
dense particles in the box. To do this, for every dark matter particle 
in the Millennium Simulation we have estimated the local density 
by using a standard (Monaghan & Lattanzio 1985; Hernquist & 
Katz 1989) smoothing kernel averaging over 32 neighbours. We 
produced a list of the most underdense particles and then radially 
sorted surrounding particles in distance from these points. In this 
way the radius of a volume underdense by a factor of 10, cen- 
tred on the most underdense particle, was calculated. All particles 
with densities less than 0.035 of the cosmic mean and further than 
2/i _1 Mpc from an even less dense particle were tried as prospec- 
tive centres. After cleaning this catalogue by removing those posi- 
tions that lay within a larger void and limiting the size to be larger 
than 5/i~ 1 Mpc in radius, this technique produced 3024 potential 
void centre candidates. We estimate that this catalogue of potential 
void centres is greater than 99% complete, given that we tried over 
200,000 random centres but only found 2 additional voids in the 
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last 50,000 (i.e. voids with central densities close to our limit of 
0.035). 

Once a list of initial positions for searching for voids in the 
simulation had been created, we used these positions to search for 
the maximum spheres that are empty of haloes with mass larger 
than a given value. To conduct this search we follow an algorithm 
based on a modification of the one presented in Patiri et al. (2006a). 
The haloes were identified using a minimal spanning tree that links 
together particles with density exceeding 900 times the background 
density (Thomas et al. 1998). This is to focus on the core proper- 
ties. In our case, we have used all haloes with masses larger than 
8.6x 1O 11 /i _1 M0 (i.e. those haloes with more than 1000 particles). 
Using this mass cut we are selecting haloes that contain galaxies 
with stellar masses similar to the ones used in the observational 
data. To characterise the final position of the void centres and their 
radii, we populate a sphere of radius R=5 h _1 Mpc, centered on 
each initial position, with 2000 random points. For every point in 
this sphere, we estimate the position of the closest 4 haloes lying 
in geometrically "independent" octants. Then we built the sphere 
defined by these 4 haloes. This is repeated for all the 2000 random 
points. As a characterization (position and radius) of the void we 
choose the biggest empty spherical region generated in the previous 
step. It is important to stress that the position of the void defined in 
this way normally does not match the position of the initial guess. 

To match the observations we select only those voids whose 
radius (as measured by the largest sphere that is empty of 
haloes greater than a given mass) is larger than 10 ft -1 Mpc. 
This cut produces 2932 voids in our box with a median radius 
of 14 h^ 1 Mpc. To explore the alignments of the haloes we 
have concentrated only on haloes with masses 8.6xl0 11 /i" 1 M Q 
<M< 8.6x 10 12 /i -1 MQ(i.e. we restrict the sample to galaxy- 
sized haloes) and located within shells beyond the surface of the 
voids. 

In addition, we have also created a subsample of haloes which 
contain a disc-dominated galaxy at their centre. These haloes are 
expected to have had a relatively 'quiet' life and have suffered 
smaller non-linear effects after turnaround such that their spin 
properties should remain matched to those of the baryonic mat- 
ter. To select these haloes we have used the semianalytic galaxy 
catalogue of Croton et al. (2006). We have created the subsample 
by selecting the brightest dominant galaxy within 200 h^ 1 kpc of 
our halo centre. We then select those haloes where the semianalytic 
code indicates that there is a galaxy brighter than Ma-<-23 mag 
that has a bulge-to-total (B/T) ratio 0<B/T<0.4 (i.e. a Milky-Way 
like disc galaxy). 

2.3 Halo shapes and spins 

To characterise the orientation of the haloes we have determined 
the orientation of their principal axes and their angular momentum 
vectors. The principal axes of the halo mass distribution are mea- 
sured by diagonalising the inertia tensor, 

N p 

lij = -j-j-^2mkXk,iXk,j (1) 

p fc=i 

where the sum is over all the particles in the halo N p , and the 
coordinates are defined with respect to the centre of mass of the 
halo of mass M/,. The resulting eigenvalues Mha 2 /5, Mhb 2 /5, and 
M h c 2 /5 are sorted by size, in descending order. The eigenvectors 
give the directions of the principal axes. 

The angular momentum vector of the halo is given by: 



L = ^2 "ifcr k x v k (2) 

k 

where the sum is again over all the particles in the halo and r k 
and Vk are the position and velocity of each particle relative to the 
centre of mass of the halo. 

The uncertainty in the position of the inertia axes and of the 
angular momentum has been evaluated by comparing these quanti- 
ties in the same haloes in two different runs that differ by a factor of 
20 in resolution. We found that we can measure the orientation of 
the angular momentum vector with an uncertainty (as provided by 
the full width half maximum of the ditribution) of ~14° when the 
number of particles in the (low resolution) halo is larger than 1000. 
The uncertainties in the inertia axes are 13° for the semimajor and 
semiminor axes and 20° for the intermediate axis. 

Once the angular momentum and the inertia vectors are evalu- 
ated, we estimate the cosine of the angle between those vectors and 
the vector joining the centre of the void with the centre of the halo 
R: 

^ cose -(^W\) (3) 



3 RESULTS 

Fig.Q]shows the probability density distribution P(cos 6) of the co- 
sine of the angles between the inertia axes (and the angular mo- 
mentum) and the vector joining the centre of the void to the centre 
of the halo. We show the results for two different shells: the shell 
located at the surface of the void with a width of 5% of the radius 
(i.e. 1 R„ oit j<R<1.05 R VO id), and a shell located well beyond the 
surface of the void at 1.2 R„ o4 d<R< 1.4 K VO id- We do this to high- 
light the effect of moving farther away from the void surface. The 
dashed line in this figure corresponds to the probability distribution 
of randomly distributed angles (i.e. P(cos 0)=l). 

To the best of our knowledge there is no theoretical predic- 
tion for the probability density distribution of the angles between 
the inertia axes and the vector joining the centre of the void to the 
centre of the halo. For this reason, we have used the following sim- 
ple analytical expression motivated by the planar symmetry of the 
problem (Betancort-Rijo & Trujillo 2007; in preparation) to quan- 
tify the strength of the signal: 

P ^ d ^ (l + /-l),)3/ 2 < 4 > 

If p=l, we obtain the null hypothesis (i.e. P(cos 8)=l). Values 
of p>l imply that the inertial axis tends to be aligned with the 
surface of the shell, whereas p< 1 implies that the axis tends to be 
perpendicular to the void surface. This simple analytical expression 
produces good fits to the observed distribution with reduced X 2 ;$l 
in most of the cases. The results of our fits are summarized in Table 

ni 

We find significant alignments of the inertia axes within the 
shells of the voids. The major axes of the dark matter haloes tends 
to lie parallel to the surface of the voids (i.e. there is an excess of 
haloes with large values of 8). For the minor axis the alignment 
is contrary to the major axis, there is an excess of haloes with the 
minor axis oriented in the radial direction of the voids. The inter- 
mediate axis also tends to lie parallel to the surface of the voids, 
although the signal is not as strong as in the case of the major axis. 
As expected, the signal declines in all the cases as the distance from 
the centre of the voids is increased, although this decline is slow: 
we still detect a weak signal of alignments of the major and minor 
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axes at distances as large as ~2.2 R vo id from the void centre. The 
distribution of angular momentum vectors, however, is compatible 
with a random orientation. 

To test the reliability of our results we have repeated our anal- 
ysis locating the centres of the voids at random positions within 
the whole volume of the simulation. As expected, we recovered for 
all the cases a signal which is compatible with the null hypothesis. 
To run this test we have used exactly the same number of random 
centres as the number of voids we have. The number of haloes and 
their distances to the centres of these 'fake' voids are similar in this 
control experiment to the real case. 

As stated in the Introduction, any signal of alignment in the 
angular momentum of the haloes is expected to be erased after the 
turnaround by non-linear effects. For this reason, if this signal is 
still present nowadays it should be found in haloes which have had 
the 'quietest' lives since turnaround. To explore this, we have re- 
peated the same analysis as before but this time using the subsam- 
ple of haloes with a disc-dominated galaxy at their centre. The re- 
sults are shown in Fig.|2]and the strength of the signal is quantified 
in Tabled 

Contrary to the result obtained using the full sample of haloes, 
the angular momentum vectors of the haloes with a disc-dominated 
galaxy at their centre tend to lie parallel to the surface of the void. 
To test the reliability of this signal we have run different statistical 
tests: the departure of the average (Avni & Bahcall 1980) of cos# 
from 0.5 (i.e. the expected value in the null hypothesis case) and 
the Kolmogorov-Smirnov (K-S) test. Both tests reject the null hy- 
pothesis at the 99.8% level. Our results are robust to changes in 
the ratio B/T ranges selected (i.e. we still find a significant signal 
for the alignment of the angular momentum including haloes with 
galaxies with B/T<0.6). In addition, we have also checked that in- 
cluding those haloes contained in our list of voids with R<10 h~ 
Mpc do not alter (within the error bars) our results. This is as ex- 
pected beacause almost all of our initial void centres produce voids 
larger than 10 h _1 Mpc in radius. 

It is worth noting that the maximum signal is found when we 
select a shell of width 1 R„ oit j<R< 1-07 K VO id- For this case, the 
null hypothesis is rejected at 99.9%. On the other hand, the inertia 
axes show the same trends in this subsample as in the previous case 
using all the haloes. 

To characterize the strength of the alignment of the angular 
momentum we have followed two approachs. We have used (as be- 
fore) Eq. (4) and, secondly, we have compared our result with the 
theoretical prediction for this quantity from Lee (2004) within the 
framework of the TTT The strength of the intrinsic galaxy align- 
ment of the galaxies with local shear at the present epoch is ex- 
pressed as the following quadratic relation (Lee & Pen 2002): 

1 + c - - 

< LiLj >— — - — 5ij — cTikTkj, (5) 

where L is the halo angular momentum (spin) vector and T is 
the rescaled traceless shear tensor T defined as Tij=Tij /|T| with 
Tij=Tij—Tr(T)8ij/3, and c is a correlation parameter introduced 
to quantify the strength of the intrinsic shear-spin alignment in the 
range of [0,1]. To estimate c we have used the analytical approxi- 
mation suggested in Lee, Kang & Jing (2005): 

P( M ) = (l-|)+|(l-M 2 ) (6) 

The values of the parameter c we obtain in the different shells 
are summarized in Table [T] When c=l the strength of the galaxy 
alignment with the large-scale distribution is maximum, whereas 



c=0 implies galaxies are oriented randomly. In the inner shell, we 
measure in this case c=0. 151 ±0.046. It is worth noting that the 
value of c measured in this work is a lower limit of the true value 
because it has been evaluated without any attempt to correct for the 
smoothing produced by the uncertainty in measuring the angular 
momentum vector of the haloes. Consequently, the strength (and 
statistical significance) of the observed alignment should be higher. 
To quantify this effect, we have re-estimated c again but this time 
using the theoretical prediction convolved with our error function 
in measuring the angular momentum. After fitting the convolved 
function to the data, we obtain c=0. 158±0.045. 



4 DISCUSSION 

In this paper, we have shown that when haloes are selected in order 
to contain a Milky- Way like disc galaxy at their centres the angular 
momentum of the dark matter halo is oriented preferentially paral- 
lel to the surface of the voids. Observationally, the same alignment 
is detected using the baryonic matter (Trujillo et al. 2006). These 
two pieces of information are in agreement with the TTT predic- 
tion that both the dark and the baryonic matter component have 
conjointly acquired their angular momentum before the moment of 
the turnaround. Interestingly, the signal detected in the real obser- 
vation c=0.7^o'2 is higher than the one found in the simulations 
c=0. 151 ±0.046. This is to be expected taking into account that the 
signal in the dark matter haloes should be erased by non-linear ef- 
fects such as exchange of angular momentum between the haloes. 
Future work, consequently, should explore the strength of the align- 
ment of the haloes at the moment of turnaround. At that early epoch 
the strength of the signal should be as strong as the one measured 
using the baryonic component. Porciani et al. (2002b) shows hints 
that this should be the case by comparing the relation between the 
halo spin and the linear shear tensor at different redshifts from z=50 
to z=0. 

We have compared our work with previous analysis of the 
alignment of the inertia axes and angular momentum with their 
surrounding large-scale structure using different simulations. We 
concentrate first on the alignment of the angular momentum: using 
the void framework, neither Heymans et al. (2006) nor PA06 have 
found a signal of the alignment of the angular momentum within 
the shell of the voids. Their results are easily understood taking 
into account that no preselection of the haloes was done in either 
of these works and that they explored the signal in just one shell 
of width 4 h^ 1 Mpc. In fact, if we mimic their analysis we find 
c=-0.014±0.012 (p=0.996±0.005), in agreement with their find- 
ings. In the same shell, selecting those haloes with 0<B/T<0.45 
produces c=0.030±0.021 (p=1.012±0.008). This is a factor of 5 
weaker than the signal we find in the closest shell to the void sur- 
face. As we have seen in this work, the signal is only significant at 
the void surface, consequently, taking a wide shell can mask it. In 
the particular case of Heymans et al. (2006) their mass resolution 
limit could be another source of concern, since it is a factor of ~10 
worse than the present simulation. 

Comparison with other work is more complicated since the 
analysis of the alignments is done using a different scheme. BS05 
measure the alignment of the angular momentum along filaments. 
Interestingly, they found that the angular momentum of galaxy 
mass haloes shows a weak tendency to point along filaments, while 
those of group and cluster mass haloes show a very strong ten- 
dency to point perpendicular to the filaments. The significance 
and strength of their signal in terms of the c parameter is, how- 
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Figure 1. Probability density distribution of the cosine of the angles 8 between the inertia axes (and angular momentum) and the vector joining the centre of 
the voids with the halo centres. The error bars on each bin is the Poissonian error and (to avoid confussion) are only plotted for the innermost shell bins. The 
null hypothesis (i.e. a sine distribution) is represented by the dashed line. 



Table 1. The strength of the alignments on the different shells. To measure the alignment of the inertia axes we have used the parameter p (p=l in the null 
case), and c (Lee 2004) in the case of the angular momentum (c=0 in the null case). 



Shell 


Maj. A. 


Int. A. 


Min. A. 


A. M. 


A.M. 


Number of 


R void Units 


P 


P 


P 


P 


c 


Haloes 


Total Sample 


1.00 <R<1.05 


1.218±0.012 


1.061±0.011 


0.778±0.008 


1.007±0.010 


0.004±0.026 


13128 


1.05 <R<1.10 


1.165±0.014 


1.060±0.013 


0.800±0.010 


0.994±0.012 


-0.014±0.032 


8644 


1.10 <R<1.20 


1.150±0.009 


1.057±0.008 


0.820±0.007 


0.992±0.008 


-0.020±0.020 


21567 


1.20 <R<1.40 


1.111±0.005 


1.040±0.005 


0.863±0.004 


0.989±0.005 


-0.033±0.012 


60908 


1.40 <R<1.80 


1.070±0.003 


1.021 ±0.003 


0.914±0.002 


0.989±0.003 


-0.027±0.007 


200201 


1.80 <R<2.60 


1.030±0.001 


l.OlliO.001 


0.959±0.001 


0.995±0.001 


-0.014±0.004 


707760 


2.60 <R<3.20 


l.OlliO.001 


1.003±0.001 


0.985±0.001 


0.998±0.001 


-0.004±0.002 


1301720 


Disc-Dominated subsample 


1.00 <R<1.05 


1.172±0.021 


1.056±0.019 


0.803±0.015 


1.062±0.019 


0.151±0.046 


4212 


1.05 <R<1.10 


1.156±0.026 


1.040±0.023 


0.812±0.019 


0.994±0.022 


0.013±0.059 


2587 


1.10 <R<1.20 


1.148±0.017 


1.035±0.015 


0.837±0.012 


1.003±0.015 
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ever, not quantified. Consequently, the agreement between ours and 
their work can be done only qualitatively. Porciani et al. (2002b) 
explored the alignment of massive haloes (M>lO 13 /i _1 M ) be- 
tween their final spin and the initial shear tensor at the halo posi- 
tion. They found c=0.07±0.04. The mass regime explored by these 
authors and their comparison between an initial and a final property 
of the haloes prevents us making a direct comparison between their 



and our measurements of the c parameter. Finally, in a recent paper, 
Aragon-Calvo et al. (2006), using galaxy mass haloes find that the 
strength of alignment of their spins in walls is c=0.13±0.02. This 
is in excellent agreement with the value reported in this paper. 

If we focus our attention on the alignment of the inertia axes, 
we find that our results are in good agreement with those of BS05, 
ACC06 and PA06. All these authors find that the tendency of 
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Figure 2. Same than in Fig.[T]but this time using only haloes which contain a disc-dominated galaxy at their centres. Note that using this halo subsample the 
angular momentum vector tends to lie parallel to the surface of the void. 



the minor axis to lie perpendicular to large-scale filaments is the 
strongest of the alignments. It is interesting to note that the strength 
of these alignments seems to be dependent on the mass of the 
haloes, being stronger for the most massive (cluster-sized) ones. 
BS05 and ACC06 suggest that the different strength could be re- 
lated to fact that most massive haloes receive a larger infall of mat- 
ter from filaments. This could also help to explain the tendency 
of the angular momentum of the most massive haloes to be per- 
pendicular to the filaments. In this sense, the cluster mass haloes 
would acquire most of their current angular momentum from ma- 
jor mergers along the filaments, whereas the angular momentum of 
the galaxy mass haloes will still have memory of the initial tidal 
fields. 

The alignment of the haloes with their local large-scale struc- 
ture is not only of interest to constrain models of galaxy formation, 
it could also be relevant to explain other observational features. For 
example, the tendency of satellite galaxies to avoid orbits that are 
coplanar with their host spiral galaxies (know as the "Holmberg ef- 
fect") found in observations (Holmberg 1969; Zaritsky et al. 1997; 
Sales & Lambas 2004, but see Agustsson & Brainerd 2006) and in 
simulations (Zentner et al. 2005; Libeskind et al. 2006). This could 
be due to the preferential accretion of satellites along filaments, that 
we have seen are preferentially aligned with the major axis of the 
host halo. 

Finally, it is worth pointing out the potential importance of 
the alignments we have discussed here to strong and weak lensing 
studies. In particular, these alignments could contribute to the weak 
lensing signal producing a shear-ellipticity correlation (Hirata & 
Seljak 2004). The degree of contamination that these alignments 
will produce in the weak lensing surveys should be explored in 
future work. 
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